Semiclassical Construction of Random Wave Functions for Confined Systems 
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We develop a statistical description of chaotic wavefunctions in closed systems obeying arbitrary 
boundary conditions by combining a semiclassical expression for the spatial two-point correlation 
function with a treatment of eigenfunctions as Gaussian random fields. Thereby we generalize 
Berry's isotropic random wave model by incorporating confinement effects through classical paths 
reflected at the boundaries. Our approach allows to explicitly calculate highly non-trivial statistics, 
such as intensity distributions, in terms of usually few short orbits, depending on the energy window 
considered. We compare with numerical quantum results for the Africa billiard and derive non- 
isotropic random wave models for other prominent confinement geometries. 
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In mesoscopic quantum systems at low temperatures 
and far away from any phase transition, many of the rel- 
evant physical phenomena can be described in the mean 
field approximation. In this scheme the excitations of the 
system are considered as a set of independent quasi par- 
ticles with energies in a small range around the Fermi en- 
ergy, which for many-particle systems is much larger than 
the single-particle ground state energies [1]. Hence, in 
this semiclassical regime, characterized by (Fermi) wave 
lengths considerably smaller than the system size, chal- 
lenges to theory are posed owing to the arising complex- 
ity of the single-particle wave functions involved. In view 
of the correspondence principle their structures depend 
sensitively on phase space properties of the correspond- 
ing classical system [2] . This has called for an increasing 
theoretical investigation of statistical properties of eigen- 
states [3] since the seminal works by Berry [4] and by 
McDonald and Kaufman [5] . This is of more than theore- 
tical interest as fluctuations of wave function amplitudes 
govern a variety of physical processes such as, e.g., pho- 
todissociation of molecules and the measured statistics of 
conductance peaks [6] in the Coulomb blockade regime 
[7]. Moreover, advances in scanning probe techniques 
and microwave experiments allow to directly uncover the 
spatial structure of waves on mesoscopic scales [12]. 

To mimic the statistical properties of wavefunctions 
in classically chaotic quantum systems, Berry conjec- 
tured [4] that chaotic wavefunctions behave as Gaussian 
random fields, and arguments coming from semiclassics 
[4], quantum ergodicity [13], and information theory [14] 
support this Gaussian hypothesis. When supplemented 
with a Bessel-type spatial two-point correlation func- 
tion, the resulting theory is known as Berry's Random 
Wave Model (RWM), since it is equivalent to consider the 
wavefunction as a random superposition of plane waves 
with locally fixed wavenumber magnitude. The RWM 
provides universal, system-independent results consistent 
with random matrix theory. It constitutes the most 
widely used statistical description of chaotic eigenfunc- 
tions, as it has been extremely successful in predicting 



bulk or spatially averaged quantities. However, obvi- 
ously, the RWM does not account for effects of confine- 
ment potentials which pose additional constraints to the 
wave functions, reducing their randomness particularly 
in the spatial region close to the boundaries. This fact 
strongly diminishes the range of applicability of the usual 
RWM, since in many experimental situations the behav- 
ior of the wave function close to the boundary is particu- 
larly relevant (e.g. when measuring tunnel rates, the local 
density of states at surfaces or boundaries, or the conduc- 
tance by attaching leads). Hence, very recently several 
papers appeared, where boundary effects have been in- 
corporated into RWM approaches, however for very spe- 
cific geometries [15-18] only or in a qualitative way [11]. 

In this Letter we construct a RWM which allows to in- 
corporate boundary effects of arbitrary confinements in- 
cluding Dirichlet-, Neumann-, and mixed boundary con- 
ditions in both billiard and smooth systems. We combine 
the Gaussian conjecture for eigenfunction statistics with 
a semiclassically exact representation of the spatial two- 
point correlation function. This enables us to account for 
confinement-induced random wave correlations in terms 
of usually few classical paths, generalizing and improving 
ideas presented in [4,19-21]. We illustrate the generality 
and strength of our technique for different systems in- 
cluding those treated in [15-18]. 

Defining the ensemble. We focus on two-dimensional 
clean, closed systems with time reversal symmetry [22]. 
We consider energy averages over a set of Nyy normal- 
ized solutions il> n (r) of the Schrodinger equation with 
non-degenerate eigenvalues E n lying inside an interval 
W = [e — Se/2,e + Se/2}. We assume Se/e <C 1, which 
can always be achieved in the semiclassical limit we are 
interested in. Considering such energy averages is stan- 
dard for disorder-free mesoscopic systems as it allows for 
random matrix approaches [23]. Moreover, experiments 
often involve averages over finite energy windows [24]. 
In particular, the energy-averaged eigenfunction inten- 
sity to be considered is proportional to the local density 
of states, relevant to many experiments such as photoab- 
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sorption, quantum transport, and ionization processes. 

At a fixed position f = (x, y) we will probe wave 
function amplitudes by means of a function F(u n ) = 
F(ip n (r)) which fluctuates when varying E n and the cor- 
responding state ip n inside W. We define the spectral av- 
erage of F at r as T{f) = Y.E n fW ^(^(0). A typ- 
ical example is the distribution of intensities, I(w:f) = 
Eb„ g w S ( w ~ I ^« (f) 1 2 ) ■ Thc definition is easily gen- 
eralized to higher-order statistics such as the spatial cor- 
relation of the intensity distribution, Y(w\,W2',fi,f2) = 
W^J2 Erl ew d ( w ^ ~ \^n{ri)\ 2 ) 5{w 2 - |^„(r 2 )| 2 ), and to 
functions F(u) depending not only on the eigenfunctions 
but also on their derivates of any order: 
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F(ui(n),..,u^(f M )) (l) 



where u*(n) = d l «.d™°"il) n {xi, j/j) with integers l a ,m a . If 
there are J different positions among the set fi , . . . , Tm > 
we call T{f\, . . . , Fm) a J-point statistics. In this paper, 
a central quantity is the two-point correlation function 
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since the average of any expression bilinear in the wave- 
function can be expressed through this correlation. 

The Gaussian conjecture. Introducing the joint prob- 
ability distribution P(u) = j^- J2e„ew — tnc 
statistics (1) can be cast into the more familiar form 
•F(ri, . . . , fju) = J_ QO F{u)P(u)du. The Gaussian con- 
jecture for the statistics of eigenfunctions of classically 
chaotic quantum systems claims that the energy ensem- 
ble is described as a Gaussian stationary process. More 
precisely, this means to assume (in the weak sense) 
P{u) = (2tt)- m / 2 (detC) _1/2 cxp(-iw(C- 1 )w), where 
the correlation matrix C = C(fi,...,fj) has entries 
c a .fi = ]^EE„gw u n u n- Since all these entries consist 
of averages over quantities bilinear in the eigenfunctions, 
the knowledge of the two-point correlation function (2) 
completely determines, under the Gaussian assumption, 
the matrix C and all statistical properties. 

Applying this approach to the intensity distribution 
I(w;r), the matrix C reduces to a single entry ci,i = 
R(f,f). Using the above expression for P(u) we find 
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Due to the presence of the boundary, R(f, f) will gener- 
ally depend on f (as will be discussed in Fig. 1) . This con- 
stitutes a non-isotropic generalization of the (isotropic) 
Porter-Thomas distribution, given by R(f, f) = const. 

The correlation of thc intensity distribution, Y, in- 
volves a 2 x 2 correlation matrix with elements Cij — 
R(fi,fj). The Gaussian integrals then give 
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Cl,lW2 + C 2 , 2 Wl 

2detC 



, (4) 



which is the non-isotropic generalization of the distribu- 
tion studied in [19,20]. 

Semiclassical construction of the correlation matrix. 
The above scheme critically depends on how precisely 
R{fi,fj) can be calculated. This is a serious issue in the 
theory of chaotic quantum systems where no analytical 
expressions for the eigenfunctions exist, and approximate 
methods are required. It turns out convenient to express 
R(fi, fj) through the Green function G(fi, fj;E + i0 + ), 

R(r i ,r j ) = --^- [ + ' lmG(r l ,r J ;E + iO+)dE, (5) 

7T Nw J e -5e/2 

since a variety of approximations exists for G. 

We start from the the exact multiple reflection expan- 
sion of the Green function [25] and consider the two lead- 
ing terms, G ~ G (0) + G (1) , to calculate R(fi,fj). 

The term G^ denotes the contribution from the direct 
path joining r*j and fj. The corresponding isotropic con- 
tribution R ls (fi, fj) to R can be calculated directly from 
Eq. (5) by means of the short-time propagator for direct 
paths. For small distances q = \f — fj\ [26] it is evalu- 
ated at the mean potential V(Q) for a local wave number 
hk = \2m{e-V{Q))] 1 /' 2 with mass to and Q = (f l +f :j )/2: 



-J (kq)T 
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Here, T(x) = smx/x is a window function; T(x) = 1 
corresponds to Berry's celebrated result [4,21,27] for the 
isotropic RWM. By choosing be — H/ti with r; the bal- 
listic time scale associated with the system size I, one 
obtains T(2q/l) [28], i.e. the correlation function is sup- 
pressed on distances of the order of the system size [29]. 

The second term, G^ , represents all quantum paths 
between r, and fj hitting the boundary once (including 
non-specular reflections) . 

The power of the representation (5) for the correlator 
R(fi,fj) combined with the Green function expansion is 
demonstrated for the Africa billiard [30] depicted in the 
left inset of Fig. 1. The numerical evaluation of R(f 1 f) 
(see Fig. 1) and R(0,f) (right inset in Fig. 1) within 
this approximation is extremely fast and the results (thin 
lines) show considerable agreement with numerically ex- 
act, but time consuming quantum mechanical reference 
calculations (symbols). The boundary effects (e.g, the os- 
cillations in R(f,f}) are adequately incorporated in the 
one-bounce treatment, but evidently beyond the range 
of applicability of the isotropic RWM, Eq. (6) (dashed 
lines), which yields, e.g., R(f,f) = const. 

In the semiclassical limit the terms in the multiple re- 
flection expansion can be further approximated by the re- 
spective semiclassical Green function [2] G sc (n , fj ; E) = 
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(ih\ / 2nih)^ 1 iDy] 1 / 2 exp (j-S 7 — ^ 7 §) where 7 now 
labels the classical paths joining f with fj. Z? 7 and 
/z 7 are smooth classical quantities, and S 1 {fi,fj] E) = 
J p- dq is the classical action along the path. For energy 
windows satisfying Se/e <C 1, the energy integral (5) then 
yields the two-point correlation function 



R(fi,fj) 



2mde ^ / T 7 

TW 



(7) 



x cos [Sj(ri, fj\ e)/h — /i 7 7r/2] 



in terms of classical paths. In Eq. (7), T 1 is the traversal 
time of path 7, and tw = 2h/Se is a characteristic cut-off 
time associated with the energy window W. Eq. (7), rep- 
resenting a generalization of the correlator R conjectured 
inRef. [21] (given by taking T(x) = 1), incorporates three 
significant advantages: First, it is semiclassically exact. 
Second, it allows to appropriately describe the statistics 
for a chosen energy window by controlling the longest 
path to be included via tw 1 while this time scale is miss- 
ing in the case T(x) = 1. Third, most importantly, it is 
compatible [31] with the definition of a correlation, (2), 
contrary to the correlation used in [21]. 




FIG. 1. Two-point correlation function R(r,f) and R(0,f) 
(right inset) for r pointing along the line indicated in the 
Africa billiard (left inset). The symbols mark numerical quan- 
tum results for 7?, Eq. (2) [32]. The thin lines depict the 
semi-quantum prediction employing Eq. (5) where the Green 
function is approximated by a sum over paths, including 
diffraction effects, with at most one reflection at the bound- 
ary. The dashed lines show the isotropic RWM result (6). 

The isotropic correlation i? ls , Eq. (6), turns out to be 
extremely robust with respect to an additional spatial 
average. To see this we note that for fixed Q the inte- 
gration over the relative position q in any small region 
will contain the continuous set of paths joining fi with 
f + q directly and the contribution from non-direct paths 
being isolated in chaotic systems. Hence in the semiclas- 
sical limit the spatial integration over the continuous set 
of direct paths yields the dominant contribution which 
coincides with the isotropic result. 



On the contrary, for a pure energy average the con- 
tribution from non-direct paths to the correlation is of 
the same semiclassical order than that from direct paths. 
However, the width 5e (corresponding to the number N\y 
used to define the ensemble) determines the maximum 
length of the non-direct paths contributing to the cor- 
relation function. The major step beyond the isotropic 
case is then to consider an energy window such that only 
the direct and shortest non-direct paths significantly con- 
tribute to the correlation function, i.e. a situation which 
is also particularly experimentally relevant. 

To this end we need to specify the non-direct paths 
more precisely. In billiard systems the first non-direct 
contribution to R is given by a sum £ p R^ (f , fj ) over 
usually few classical trajectories p hitting the boundary 
once. For given initial and final positions r*j, fj each one- 
bounce path p is uniquely characterized by the position 
f p where it is reflected. The path length is L p = L ip + L 



with Lj 



f p \, Lj p 



j — f p \ . Denoting by k p and 
6 p the local boundary curvature and reflection angle at 
f p , a simple calculation yields for each path 



R^(f,f 3 )^ri kL ^ e 



2k Tj 



LipLjp 
LpCosOp 



- 1 



Ay/2TrkL p 



cos ( kL v 



(8) 



to the correlation function. Here A is the billiard area, 
and (j) p takes into account the boundary conditions at the 
reflection point, as given e.g. in [33]. 

The function R(n,fj) = R is {fi,fj) + £ p i? (p) (^, ^0, 
together with the semiclassical expressions (6) and (8), 
provide the entries for the correlation matrix C, from 
which arbitrary statistical measures (such as /, Eq. (3), 
and T, Eq. (4)) for the wave functions can be deduced 
beyond the isotropic case. Moreover this semiclassical 
correlation yields closed analytical expressions for sta- 
tistical quantities for chaotic systems as the stadium-, 
cardiod-, or Sinai-billiard, since in these cases all the pa- 
rameters required are readily calculated from geometrical 
considerations. 

The preceding discussion is easily generalized to sys- 
tems with Aharonov-Bohm flux lines or with smooth 
boundary potentials, where the first non-direct contri- 
bution includes paths with one classical turning point. 

Non-isotropic random wave models: one path is 
enough. In the following we demonstrate the power of the 
general semiclassical scheme outlined above by comput- 
ing wave function correlators for selected, representative 
examples. First we show the role of confinement effects 
by considering points f close to the boundary of a bil- 
liard system. In [15] this situation is treated by approx- 
imating the boundary by an infinite straight line y = yo 
and considering an ensemble of random superpositions of 
plane waves ip r (f) satisfying the general, mixed boundary 
condition (d y ip r (f) sin a + kip r (f) cos a)\ y=ya = 0. Here 
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a is a generally position-dependent parameter and k is 
the local wavenumber. By ensemble average a variety 
of two-point correlations was derived in [15] and used 
to calculate specific statistical observables. To illustrate 
our method and for the sake of comparison we consider 
in detail just one such average, namely (ip r (r)d y ip r (f)) . 
In terms of the two-point correlation function it reads 
(il>(r)d y ip(r)) = ^(dy.+dy^Rifufj^y^y^y. Closcto 
the boundary only the direct and the shortest non-direct 
path contribute. For mixed boundary conditions the ex- 
tra phase <j)p in Eq. (8) is given semiclassically [33] by 
(j) p — 7r — 2arctan(tanacos6> p ). Substitution into Eq. 
(8) gives the approximate correlation function close to 
an arbitrary boundary. To leading order in k we get 



<V>(W(r)> 



2kd(r)Se 



1 



sin 



(2kd(r) - 2a - ^) , 



(9) 



where d(r) is the shortest distance from f to the bound- 
ary. In the limit of flat boundaries or very short distances 
d{r) the semiclassical results represent the large-fc limits 
of the integral expressions given in [15]. The Dirichlct 
and Neumann cases considered in [16] are particular cases 
of Eq. (9) corresponding to a = and a = tt/2. 

To show how to apply the semiclassical scheme for 
more general situations, we consider now a smooth po- 
tential barrier. In [17] an ensemble of random superpo- 
sitions of Airy functions in y-direction and plane waves 
in x-direction is introduced, which satisfies locally the 
Schrodinger equation for the potential V(r) = Vy. En- 
semble average then gives, up to an overall constant, 
(r(r)d v r(r)} = /~Ai[*(j/,Q)]Ai'[*(y,Q)]dQ with 

= {Vh- 2 Y (y - yo) + (Vft- 2 )~*Q 2 , where 
ya = e/V is the classical turning point. Ai(x) and Ai (x) 
is the Airy function and its derivative [15,17]. The clas- 
sical paths (with no or one turning point) required to 
construct the corresponding average via the semiclassical 
correlation function can be calculated in closed form as 
they are just parabolic flights. Using Eq. (7) and keeping 
only terms to leading order in ft -1 , we finally get 
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(2mV{y -y)) 1/2 Se 



nv 



{2irhf/ 2 N w 
1 . / Ay/2mV . 



sm 



2hVW(y - y) 2 \ 



3/2 



which is again the asymptotic limit of the integral ex- 
pression presented above. Correspondingly, we recover 
the asymptotic limits of the results for the geometries 
studied in [18] in terms of a small number of paths. 

To summarize, we showed how to efficiently treat wave- 
function statistics for closed systems by merging statis- 
tical with semiclassical concepts. We demonstrated that 



all known (to us) results [15-18] for specific, non-isotropc 
Random Wave Models are particular cases of the general 
approach presented here. It provides closed analytical ex- 
pressions for statistical measures in terms of geometrical 
quantities and builds the framework for incorporating ar- 
bitrary boundary conditions and confinement geometries. 
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